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Abstract 

We present a simple method for incorporating the surface tension effect into an iterative 
conformal mapping model of two-dimensional diffusion-limited aggregation. A curvature- 
dependent growth probability is introduced and the curvature is given by utilizing the 
branch points of a conformal map. The resulting cluster exhibits a crossover from compact 
to fractal growth. In the fractal growth regime, it is confirmed, by the conformal map 
technique, that the fractal dimension of its area and perimeter length coincide. 
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1 INTRODUCTION 

Since its introduction by Witten and Sander in 1981^, the diffusion-limited aggregation (DLA) 
has been one of the most important models of growth processes which are dominated by diffusion 
process or, in a wider sense, of non-equilibrium and non-linear statistical physics and stochastic 
processes. The dynamical rule of this model is very simple: Initially a seed particle is placed 
at the center. A random walker is released far from the cluster and when once it touches the 
cluster, it stops walking and becomes a part of the cluster. Then next walker is released, 
The process repeats again and again. 

The cluster grown under the above-described process forms a fractal structure having sev- 
eral, typically five for a two-dimensional cluster, main branches and each branch has many 
subbranches hierarchically. Growth patterns like this have been observed for many kinds 
of experiments: electrodepositiorP, viscous fingering^ ^1^, crystal growth with its anisotropy 
suppressed^, bacteria colony growthP, and so forth. These experiments, computer simulations, 
and theoretical studies have clearified the properties of the DLA, but some fundamental prob- 
lems to be understood remain. 

About fifteen years ago, a novel approach was proposecP for two-dimensional DLA. This 
approach utilizes iterative conformal map, which automatically satisfies the Laplace equation, 
since the DLA is dominated by the field which satisfies it. This approach has led to many 
interesting results: For example, scaling laws associated with the Laurent expansion of the 
mapP, higher precision multifractal analysis and the role of randomness^. 

It is well known that a flat interface becomes unstable in a diffusion-dominated growth 
process^ The instability is suppressed within a certain length scale by the surface tension 
effect which stabilizes the flat interface l4 . Therefore the incorporation of this effect is important 



2 



for considering the validity of a model being compared with experimental results and the relation 
between the randomness of a growth rule and regular growth patterns. There have been some 
attempts to incorporate the effect into the original form of the DLA, by introducing a sticking 
probability dependent on the curvature of the growth point Oil ED anc j m0 re ingeniously, by 
reconstructing the interface by random walks between two points on the surface^. 

In the present article we present how to incorporate the surface tension effect into an iterative 
conformal mapping approach. We follow the idea of a probability dependent on the local 
curvature at the growth point?™. The curvature is calculated by utilizing the properties of 
the conformal map. 

The rest of the article is organized as follows: In Sec. 2, we describe the model. We briefly 
review the iterative conformal mapping approach. And then we introduce a probability with 
curvature dependence and present how to calculate the curvature. The patterns formed by this 
model are shown and qualitatively discussed in Sec. 3. A scaling analysis is given in Sec. 4. 
Sec. 5 is dedicated to the summary. 

2 MODEL 

2.1 Iterative Conformal Mapping Model of DLA 

Let & n >(z) be the conformal map which maps the external region of the unit circle \z\ = 1 in 
the complex plane C onto the external region of the two-dimensional n— particle cluster in the 
physical space. The image of the unit circle under & n >(z) is the interface of the cluster. 

First let (f>g t \ be the elementary map which makes a semicircular "bump" of size at 
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Then $( n ) (z) is constructed by iteration of the elementary maps: 

$ (n) W = ^o^°-o^(4 (3) 

Note the order of the composition. 

The angle at which the n-th bump is generated, 9 n , is an independent stochastic variable 
obeying the uniform distribution on [0, 2-rr]. The bump size parameter A n is chosen as 

A" = TjfJ I it) M 2 ' (^) 

l $ n-i( e ™)r 

in order to fix the bump size in the physical plane, to the first order, where Ao is a typical bump 
area and is specified in advance. 

In FigJT]a typical cluster of 10000 bumps and Ao = 0.1 constructed by the above approach 
is shown. 

2.2 Surface tension effect 

In order to incorporate the surface tension effect into the conformal mapping model, we in- 
troduce the growth probability P(z) that the next ((n + l)-th) bump is generated at around 
z = $( n )(e ^6, ). The probability P(z) depends on the curvature at z and is written as 

P(z) = -Ak(z) + B, (5) 

where A and B are positive constants and k(z) is the curvature at z. Since P(z) is defined as 
a probability, if P(z) > 1 we set P(z) = 1. Also, if P(z) < e we set P(z) = e where e is a 



small positive constant, in order to keep the process going on and save the calculation time. It 
is necessary for e to be smaller for larger A in order to suppress the effect of randomness which 
allows more growth at a point with larger curvature. 

The introduction of the "sticking" probability, Eq.(JS]), that a random walker released far 
from the cluster sticks to a given point, to incorporate the surface tension effect, was presented 
in on-lattice^ and off-lattice^ simulations for diffusion-limited growth process. Subsequently, 
it was showri^that the evaluation of the surface tension effect could be better by reconstructing 
the interface by considering a random walk between two points on the interface. This method 
was used in simulations of viscous fingering on a lattice^"^. It was pointed outP^that there is 
a correspondence between these two methods when we consider the short-time behavior of the 
interface not so far from flatness. 

In order to obtain the curvature at a point, we utilize the fact that the conformal map 
has branch points associated with the generations of bumps^. Each branch point has a 
preimage exp[i/3^^] on the unit circle. The argument of the preimage is characterized by three 
integer indices, k, n, and j(n). The subscript k is the step when the branch point was generated 
(i.e. the fc-th bump was generated). The subscript n is the step at which the analysis is being 
done. The superscript j(n) is the order of the branch point along the arc length, emphasizing 
that it is a function of the step n. Consider the list of the arguments of the preimages at the 
n-th step, {fll^}. For each step, the preimages move on the unit circle in the mathematical 
space keeping their physical positions fixed. When bumps overlap, some branch points are 
covered and their arguments are omitted from the list. Assume that the list at (n — l)-th step 
{ftk'n-l^} 1S available. At the n-th step, the elementary map associated with the 72,-th bump 
(j)e n ,\ n has two branch points at 

exp[ia±] = (f)0 niXn (exp[if3± n }), (6) 
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where 

Pn,n = ® n ± tal1 

in order to keep the asymptotic form (f>e n ,\ n {z) ~ z for z — > oo. On the unit circle, the domain 
[a~, a+] is covered by the bump. If an element of the list is in this domain, it does not join the 
list at the next step. Therefore the list at the n-th step „ } consists of the new elements 
/9^ n and the elements of the list at the previous step, which are not omitted by the new bump 
and are updated as 

exp^^^^Jexpf^^]), (8) 

where 4>b^\ denotes the inverse function of (j)e n ,x n - Note that the order index j'(ri) is not simply 
related to j(n — 1). Each branch point Zj at the n-th step is given as 

Zj = *W(exp[i#W]) 

= $ (fc W/4 n) ]), (9) 

where exp[i/3j( " ] and exp[z/3^] are related as 

ex P [^l ( n n) ] = 4>o! M o • • • o < ljAfe+1 (ex P [^]). (10) 

Numerically the second-line representation of Eq.(j9]) gives a much more precise result than the 
first. We can draw the outline of the cluster by connecting the images of the branch points. 

Assume that X = $( n )(e ie ) is a candidate position of the n + 1-th bump which is between 
the j- and (j + l)-th branch points. Let A, B, and C be the midpoints between the (j — 1)- 
and j-th, j- and (j + l)-th, and (j + 1)- and (j + 2)-th branch points, respectively. Then 
the curvature radius R is given as the radius of the circumcircle of AABC (see Figj2]). Let 
X = /LABC be the angle pointing to the external region. If % > 7r, the curvature is positive 
and if % < 7r it is negative. The validity of this definition of the curvature will be examined in 
the next section. 

6 



2v A n 
1- A, 



(7) 



For this model, it is known that sometimes a large flat bump is generated which seals a 
"fjord" of the cluster-. In order to prevent a flat bump from being generated, we restrict the 
number of branch points, N c , covered by a bump. However, too small N c makes it impossible for 
a bump to be generated at a point of negative curvature, since branch points are concentrated 
around such point. Hereinafter, we set N c = 6. 

3 RESULTS 

Figs. Eta)-(d) show typical patterns for several values of A. Values B = 0.5, A = 0.1 and 
iV = 10000 are common to all the calculations. The dependence on B is very weak. The 
appropriate value of e depends on A, since it should be chosen in order not to affect the results. 

The parameter A corresponds to the strength of the surface tension effect. Therefore, with 
increasing A, the structure of the cluster becomes more coarse and compact and the cluster has 
fatter branches. 

We can also observe the above-described crossover of the pattern by changing the size of 
the cluster while keeping A fixed, since the surface tension is effective within a certain finite 
length scale. This is demonstrated in FigJU where A = 2.0, B = 0.5, e = 0.0001, and interfaces 
are plotted every 750 steps. We can observe that an initially two-dimensional compact cluster 
becomes fractal through tip-splitting. 

In order for our definition of curvature to be valid, the curvature radius R should be suf- 
ficiently larger than the distance between branch points. Fig. 5 shows the distribution of the 
curvature of the point at which a bump is generated for several values of A until iV = 10000. A 
dependence of the distribution on step n is not observed. It is easily confirmed that the greater 
A becomes, the more suppressed the growth at a point with large curvature is. The character- 
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istic distance between branch points is evaluated by the typical size of bump, i/Ao ~ 0.3. Thus 
it is necessary that a growth at a point with curvature larger than ~ 3 should not be allowed. 
From FigJSJ it is concluded that for A > 0.25, our definition of curvature is valid. 

4 SCALING ANALYSIS 

In this model, some quantities can be evaluated from the conformal map. The area of the 
cluster Sn is evaluated by the number of the bumps (steps) N, since the size of the bump is 
kept fixed to the first order (see Eq.(j3])). The perimeter length Ln of the cluster is evaluated 
by the number of the branch points iV&. From the Laurent expansion of the map <f>( n >(z): 

oo 

= F^z + 4 0) + F { N j) z- j , (11) 

3=1 

the characteristic radius of the cluster of N bumps Rn is given by the first order coefficient, 
which is called the Laplace radius or conformal radius and written in terms of the bump 
size parameter {Xj} as 

TV 

^~4 i) ~ri( i+A i) i/2 - ( i2 ) 

3=1 

Note that Fjp can be taken to be positive without loss of generality. 

We expect that the scaling relations shown below hold for these quantities: 

S N ~R%°, (13) 
L N ~R% L , (14) 

where D s and Dl are scaling exponents. 

Figsj6] (a) and (b) show the plot of R^ vs. Sn and R^ vs. L N , respectively for A = 1.0, 
averaged over 20 samples. For R N > R cN ~ 5, the growth is fractal, with the exponents 
Ds = Dl 1.7. In this regime Ds and Dl are the fractal dimension of the area and perimeter 
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length, respectively. Our result Ds = Dl ~ 1-7 agrees with the result for the orig inal DLA^. 
For i?Ar < RcN i the growth is compact, = 2 and Dl = 1 (see Figj4]). If a larger A is used 
and an average over more samples are taken, R c n becomes larger and the compact growth is 
expected to be observed more clearly. 

5 SUMMARY 

We presented a method for incorporating the surface tension effect into an iterative conformal 
mapping model of the DLA. In our method, a generation of a bump depends on the curvature 
of the candidate point where this curvature is obtained by utilizing the branch points of a 
conformal map. Scaling analysis of the cluster grown under the proposed model is easily 
conducted using the conformal map technique. The grown pattern exhibited a crossover at a 
certain length scale from compact growth to fractal growth, due to the surface tension effect. 
In the fractal regime, the fractal dimensions of the area and the perimeter length coincide and 
their values agreed with those of previous results. 
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Figure captions: 

Fig.l: Typical DLA cluster constructed by the iterated conformal mapping procedure, iV = 
10000 and A = 0.1. 

Fig.2: The approximate curvature at X, a candidate of the next growth point. The thick 
grey curve is the interface, which is the image of the unit circle under <I>( n ). The branch points 
of are labelled along the arc length, (j — l),...,(j + 2), A, B and C are the 

midpoints between the (j — 1)— and j—th, j and (j + 1)— and (j + 1)— and (j + 2)— th branch 
points, respectively. The radius of the circumcircle of AABC, denoted by R, is the approximate 
curvature radius at X. In this figure, the angle x is t ne exterior angle of ZABC and so the 
curvature at X is positive. 

Fig.3: Typical clusters for (a) A = 0.5, e = 0.01; (b) A = 1.0, e = 0.001; (c) A = 2.0, 
e = 0.0001; and (d) A = 4.0, e = 0.0001. N = 10000 and A = 0.1 are common values. 

Fig. 4: Interfaces of the cluster of Fig.3(c). The contours are plotted every 750 steps. 

Fig. 5: Distribution of the curvature of growth points for several values of the strength of 
surface tension effect A. The values of e are the same as those for Fig 01 with the addition of 
e = 0.01 for A = 0.25. 

Fig. 6: (a)The log-log plot of the cluster area versus radius R]y. (b)The log-log plot of 
the perimeter length L N versus for A = 1.0, e = 0.001. In each plot the broken and dotted 
lines are guides for the eyes. 
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Figure 1: 
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Figure 2: 
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Figure 4: 
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